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■ Abstract. We study the problem of stabilization for the acoustic system with a spatially 
£S| . distributed damping. Imposing various hypotheses on the structural properties of the damping 

term, we identify cither exponential or polynomial decay of solutions with growing time. Expo- 
p I ' nential decay rate is shown by means of a time domain approach, reducing the problem to an 

■ observability inequality to be verified for solutions of the associated conservative problem. In 
addition, we show a polynomial stabilization result, where the proof uses a frequency domain 
method and combines a contradiction argument with the multiplier technique to carry out a 
special analysis for the resolvent. 
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t^j- 1 1 Introduction 

od_ 

■ We consider the following system of equations: 

t-H ■ 

u t + Vr + a u = 0, in tt x R + , 
r t + divu = 0, in x M+, 

u- n = 0, on T x R+ ^ ' ' 

u(0,x) = vP(x), r(0,x) = r°(x), x £ J7, 

where Q, is a bounded domain in M. d , d = 2, 3, with a smooth boundary T, div = V- is 
the divergence operator and a £ C°°(f2), with a > on Q and such that 

3 a_ > such that a > a_ on uj. (1-2) 

Here to ^ stands for the subset of Q, on which the feedback is active. As usual n 
denotes the unit outward normal vector along T. 
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The system of equations (jl.ip is a linearization of the acoustic equation governing 
the propagation of acoustic waves in a compressible medium, see Lighthill |1C H [T2]. 
where au represents a damping term of Brinkman type. This kind of damping arises also 
in the process of homogenization (see Allaire [1]), and is frequently used as a suitable 
penalization in fluid mechanics models, see Angot, Bruneau, and Fabrie [3J. Our main 
goal is to find sufficient condition on the initial data and the function a so that the 
solution of (jl.ip stabilizes to zero for t — > oo and, if this occurs, we are interested in the 
rate of this decay. 

Let L 2 (Q) denote the standard Hilbert space of square integrable functions in fi. To 
avoid abuse of notation, we shall write || • || for the L 2 (0)-norm or the L 2 (f2) rf -norm. 
Denoting H = (L 2 (Q,)) d x L 2 (f2), we introduce the operator 

A = ( ^ J J : V(A) = {(u, r) G H, (Vr, divu) G H, u.n\ T = 0} C H H, 



and 



B=(^\e C((L\n)) d , H),B* = {^ ) G C(H, (L\n)y 
Accordingly, the problem (jl.ip can be recast in an abstract form: 



Z t (t) + AZ(t) + BB*Z(t) = 0, t > 0, 
Z(0) = z°, 



(1.3) 



where Z = (u,r), or, equivalently, 



Z t (t) = A d Z(t), t > 0, 
Z(0) = z°, 



(1.4) 



with = -A - BB* with V{A d ) = V{A). 

It can be shown (see Section [5] below) that for any initial data (-5°,r) G D(A) the 
problem admits a unique solution 

(u,r) G C([0,oo);P(^4)) nC 1 ([0,oo);ff). 

Moreover, the solution ({?, r) satisfies, the energy identity 

E(0) - E(t) = \\V^u(s)\\ 2 mn))d ds for alii > (1.5) 



with 

FJt) = 

where we have denoted 



E{t) = \ \\(u(t),r(t))\\ 2 H ,Vt>0, (1.6) 
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Using (jl.6j) and the standard density arguments, we can extend the solution operator 
for the data (it°,r) G H. Consequently, we associate to the problem (jl.ip (or to the 
abstract Cauchy problems (jl.3j) . (|1.4|) ) a solution (semi)-group that is globally bounded 
in fl". 

As the energy E is nonincreasing along trajectories, we want to determine the set 
of initial data (vP,r°) for which 

E(t) -> as t -> oo. (1.7) 

Such a question is of course intimately related to the structural properties of the function 
a, notably to the geometry of the set to on which the damping is effective. 

In this paper, we characterize the set of initial data for which (|1.7p holds in terms of 
the set u, and, eventually we obtain some information on the rate of decay. In particular, 
we establish an observability inequality for the associated conservative system yielding 
exponential decay and use a frequency domain method, combined with the multiplier 
technique, to obtain polynomial rate of decay. It is worth-observing that the associated 
conservative system coincides with the standard linear wave equation, supplemented 
with the Neumann boundary conditions, where the asymptotic behavior of solutions is 
relatively well understood. 

The paper is organized as follows. Section [2] summarizes some well known facts 
concerning the acoustic system (jl.ip . In section [3l we examine the spectral properties 
of the generator Ad and establish a strong stability results. Section d] addresses the 
exponential and not exponential stability results. In Section [5J we prove exponential 
stability for a modified system, with a slightly different damping law, by using an 
observability strategy. Polynomial stability of the modified system is studied in Section 

m 



2 Preliminaries 

We start with a simple observation that the problem (jl.ip can be viewed as a bounded 
(in H) perturbation of the conservative system 

f u t + Vr = 0, inOxR+, ) 

\ r t + divu = 0, in O x K+, J 1 ' ' 

which can be recast as the standard wave equation 

r t ,t - Ar = 0. 

Consequently, the basic existence theory for (jl.ip derives from that of (|2.8p . 

Moreover, since the boundary T as well as the damping coefficient a are smooth, 
solutions of (jl.ip remain smooth as soon as we take it , r° smooth and satisfying relevant 
compatibility conditions as the case may be. In what follows, we may therefore deal 
with smooth solution, whereas the results for data in H can be obtained by means of 
density arguments. 
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2.1 Long-time behavior 

The operator Ad possesses a non-trivial kernel that is left invariant by the evolution, 
namely, solutions of the "stationary" problem 

Vr + au = 0, V ■ u = 0, u ■ n\ r = 0. (2.9) 

Thus the stationary field u is solenoidal, and, integrating (|2.9p over f2 yields 

{? = in supp a, Vr = 0. 

Accordingly, we introduce the space 

E = Ker[^4rf] = {(u, r) | V • u = 0, "u| S uppa = 0, u ■ n|r = 0, r = const}, 

together with its orthogonal complement (in H) denoted Hq. 
It is easy to check that 

{Ad(w, s), (u, r)) H = for any (w, s) £ V(A), (u, r) <E E; 

in particular, the solution operator associated to (jl.ip leaves both E and Hq invariant. 
Consequently, the decay property (jl.7p may hold only for the initial data emenating 
from the set -Ho- 

3 Strong stability 

The following observation can be shown by a simple density argument: 
Lemma 3.1. The solution (u,r) of (ll.lh with initial datum in D(Ad) satisfies 

E'(t) = - a\u\ 2 dx < 0. (3.10) 
Jn 

Therefore the energy is non-increasing and (jl.5p holds for all initial datum in H. 

As already shown in the previous section, the strong stability result (jl.7p may hold 
only if we take the initial data 

(u°,r°) eH = Kev[A d } ± . 

There are several ways how to show (jl.7p , here we make use of the following result due 
to Arendt and Batty [3]: 

Theorem 3.2. Let (T(t))t>o be a bounded Co-semigroup on a reflexive space X. Denote 
by A the generator of (T(t)) and by o~(A) the spectrum of A. If a (A) n iM. is countable 
and no eigenvalue of A lies on the imaginary axis, then lim T(t)x = for all x £ X. 

In view of this theorem we need to identify the spectrum of Ad lying on the imaginary 
axis. 

First we look at the point spectrum. 
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Lemma 3.3. Suppose that \u\ > 0. If X is a non-zero real number, then iX is not an 
eigenvalue of Ad- 

Proof. Suppose that 

iXu + Vr + a u = 0, 



(3.11 

iXr + div u = U. 



Prom (|3.1ip we deduce that 



iX / — \r\ )dx + / a\u\ dx = 0; 
In Jn 



whence 

^Isuppa — 0, and, consequently, u, r solve (|3.1ip with ct — 0. In particular, we 

get 

-Ar = A 2 r, Vr • n\ r = 0, r| suppQ = 0, 
and, by unique continuation for elliptic problems, we get r = 0. □ 



In accordance with Lemma 13.31 and the discussion in the previous section, A = is 
the only possibly eigenvalue of Ad on the imaginary axis. 

Next, we show that Ad has no continuous spectrum on the imaginary axis, except 
eventually zero. 

Lemma 3.4. Suppose that > 0. If X is a non-zero real number, then iX belongs to 
the resolvent set p{Ad) of Ad- 

Proof. In view of Lemma 13.31 it is enough to show that iXI — Ad is surjective. 
Hence given a vector \f,p) S H, we look for (u, r) E V(Ad) such that 



(iXI-A d )(u,r) = (f,p). (3.12) 
By the definition of Ad, we obtain 

J iXu + Vr + au = f, 
\ iXr + div u = p. 

Assuming that u and r exist we can write 

1 



(-Vr + /). (3.13) 
iX + a 

Inserting this expression in the second identity we obtain the differential equation in r: 

iXr — div f— X7r)=p — div (— f] in Q, Vr • n|r = 0. (3.14) 

yzA + a / \iX + a J 

Multiplying this identity by a test function s G ff 1 (il), integrating in and using 
formal integration by parts we get the problem: 

[ (iXrs + — - — Vr -Vs) dx = F(s),Vs £ H\n), (3.15) 
Jn\ iX + a J 
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where 

F(s)= I (ps + —^ — f-Vs] dx. 

We use the Fredholm alternative by splitting the left-hand side of (|3.15p into its 
principal part 

[I 

dn(r, s) = — Vr • Vsdx, 

pK J n iX + a 

and its lower order term 

ao(r, s) = iX / rsdx. 



The principal part is a continuous sesquilinear coercive form on 

H^ n (n) = {v€H 1 (n), [ vdx = 0}; 

hence it induces an isomorphism A p from H^fl) into (H^Q))'. The mapping Aq from 
H-^(Q) into (H^ip.))' induces by ao is clearly given by iXI and is therefore compact. 
Consequently by the Fredholm alternative, A p + Aq is an isomorphism from iJ^(O) into 
(H^iQ))' if and only if it is injective. But the injectivity of A p + Aq is equivalent to 
the injectivity of iXI — Ad- Indeed let r G H^Cl) be such that (A p + Aq)t = or 
equivalently such that 

/ (iXrs + — - — Vr • Vs | dx = 0,Vs G #1(0). (3.16) 
Jn V *A + a / 

Since r is of mean zero, this identity remains valid for all s G i? 1 (0). Hence by taking 
s G 2?($7) we find that 

iAr-div ( — - — Vr ) = in fi. (3.17) 

By taking s = r in (|3.16p . we find 

/ (iA|r| 2 + — - — |Vr| 2 ) dx = 0. 
Jn\ iX + a J 

By taking the real part of this identity, we find that 

a\Vr\ 2 = in Q. 

This implies that 

Vr = in uj. 

Hence by ()3.17p . r = in uj and by the Holmgren uniqueness theorem we deduce that 
r = 0. 

In conclusion A p + Ao is an isomorphism from H^Q) into (H^Q))' which implies 
that there exists a unique solution r G H^iQ) of 

/ (iXrs + — - — Vr • Vs] dx = F(s),\/s G i?,L(0), 
V *A + a y 

As .F(l) = 0, we deduce that r G H^(fi) is solution of (I3.15D and relation (I3.12p 
holds. □ 
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These Lemmas and Theorem 13.21 leads to 



Corollary 3.5. Let (u,r) be the unique semi-group solution of the problem lil.l]) ema- 
nating from the initial data (u°,r°) £ H. Let Pe be the orthogonal projection onto the 
space E = Ker[^] in H , and let 

(w,s)=P E (H°,r°). 

Then 

|| (it, r)(t, •) — (w, s)\\h — > as t — >■ oo. 



4 Exponential stability 



Now, we may ask, under some conditions on the damping coefficient, the convergence 
in Corollary 13.51 is exponential in time. 

Theorem 4.1. If the damping coefficient a is not uniformly positive definite, meaning 

inf a(x) = 0, 

then the system (jl.ip is not exponentially stable. Conversely, if a is uniformly positive 
definite, meaning 

oj = ft, 

then the system (ll.lh is exponentially stable, specifically, 

\\{u,r)(t, OHff < exp(-Li) || (it , r°)||#, L > 0, whenever r° dx = 0. 

Proof. According to [9] (see also [2]), the exponential stability of the system (jl.ip is 
equivalent that the undamped system, i.e. 



(f) t + Vp = 0, in ft x R+, 
pt + <iti;^> = 0, in $7 x R + , 
<j>- n = 0, on T x R+, 
0(0, x) = 0°(x), p(0, x) =p°(x), x G ft, 

satisfies the following inequality: There exist positive real numbers T, C such that 



(4.1) 







2 


0V) 


/ / ol{x) 




dxdt > C 


lo Jn 









(L 2 (n) d xL 2 (n) 



(4.2) 



V(0 t) ,p°) E (L 2 (ft)) d x L 2 (ft) such that f p° dx = 0. 

Since this estimate is well-defined in the energy space, it holds if and only if it holds for 
strong solutions. 

As the conservative system (|4.2p admits solutions that are constant in time, namely, 

u G L 2 (ft), V-n| r = 0, p = 0, 
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it is clear that exponential stability cannot hold if inf^gQ a(x) = 0. Indeed, as a is 
smooth, we can always find for any 5 > small a solenoidal compactly supported 
function cf) in fi such that 

Consequently, relation (j4.2j) cannot holds uniformly for any choice of the data. 

On the other hand, suppose that a is bounded below away from zero on the whole 
set fi. Writing the vector field u as its Helmholtz decomposition 

u = H[u] + Vtp, (4.3) 

where H denotes the standard Helmholtz projection onto the space of solenoidal func- 
tions, it is enough to verify the observability criterion (|4,2p for (j) = V<£>. In such a case, 
however, the conservative system (|4.ip reduces to the standard wave equation and (|4.2p 
is obviously satisfied as the damping acts uniformly on the whole domain fi (see also 
Section [5]). □ 



5 Changing the damping law 



As we have seen before, system (jl.ip is exponentially stable if and only if a is uniformly 
positive definite. We will show in this section that if we change the feedback law in 
order to filter the divergence free vector fields, then we will get exponential stability for 
a quite large set of a. 

In view of the Helmholtz decomposition (|4.3p . denote by P the orthogonal projection 
on the closed subspace of L 2 (fi) d 

V := lv<p : (p e H\Sl) J ip dx = oj . 

Then in (jl.ip we change the damping term a u by P(a Pu) and consider the system 



u t + Vr + P(a Pu) = 0, in x R+, 
r t + divu = 0, infix R + , 
u.n = 0, on r x R+, 
u(0,x) = v°(x), r(0,x) = r°(x), x G fi. 

Accordingly, we arrive at the system 

' Vipt + Vr + P{aVip) = 0, in fi x R+, 
r t + Aip = 0, in fi x R+, 
Vip ■ n = 0, on T x R+, 
ip(0,x) = tp°(x), r(0, x) = r°(x), x € fi, 

where 

/ ip° &x= I r° dx = 0. 
Jn Jn 
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By virtue of [9] (see also [2] ) , the exponential stability of the system (|5.2p onVxV 
is equivalent to the following property of solutions to the undamped system: 

There exist a positive real number To such that for all T > To, there exists C > 
such that 



T 



a(x) \V^(x,t)\ 2 dxdt > C ||(VV°,p )||( LWxL2(f7) (5.3) 



(5.4) 



'0 JQ 
for any ip, p satisfying 

VV>t + Vp = 0, in n x R+, 
p t + A-0 = 0, in O x R+, 

• n = 0, on T x R+, 
ip(0, x) = ip (x), p(0, x) = p°(x), i£fi, 

G H\fl), p° G L 2 (n), [ -0° dx = [ p°dx = 0. 

However, the system (|5.4|) can be written in the form of a standard wave equation 
with the Neumann Laplacian: 

iptt + A0 = 0, in O x R+, 

Vip-n = 0, on r x R+, (5.5) 

ip(0,x) = ip°(x), ipt{0,x) = if) l (x) = -p°(x), xen, 

whereas the observability inequality (|5.3p reduces to 

£ J a(x) \V4>(x,t)\ 2 dxdt> C ||(V<,^)||( LWxL2(n) . (5.6) 

Thus we have obtained the following exponential stability result: 
Theorem 5.1. The system (|5.ip is exponentially stable, meaning, 

\\(u,r)(t,-)\\ H < exp(-Lt)||(tt°,r )||H for a certain L > 0, 

whenever 

if = V(p°,p° € H\Q,), [ <p°dx = [ r°dx = 0, 



if and only if any solution ip of the wave equation (|5.5p satisfies the observability in- 
equality (|5.6p . 



Validity of the observability inequality (|5.6p is related to the general discussion in 
Bardos et al. [5] (see also Zuazua [13 [18]). Here we note that (j5.6[) holds provided 
u) contains a neighborhood of the whole boundary T, in the sense that there exists a 
neighborhood O of T in R rf such that n O C uj. First we notice that by the arguments 
of section 4 of |15| with a = and M{u) = mVu + ^-^-u (see estimate (4.10)), where, 
as usual, m(x) = x — xq for some xq G R rf , there exits C > such that for all T > 0, we 
have 

2T£(0) < C(E(0) + f [ {\^ t \ 2 + \V^\ 2 + \^\ 2 )dxdt, (5-7) 
io Jnno' 
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where O' is a sufficiently small neighbourhood of the boundary (such that O'cO). On 
the other hand, fixing a cut-off function r\ such that t] = 1 on O' and with a support O, 
we can write 

|V^| 2 dx(i£ < / / r/\ip t \ 2 dxdt. 
'o ^nno' ■/ o Jn 

Hence by integration by parts in time, we get 

/ \ip t \ 2 dxdt < — I I rjiptt'ipdxdt + / r]ip t 'ipdx 
o Jnno' Jo Jn Jn 



Now using the first identity in (|5.5p and an integration by parts in space, we get 



\ipt\ 2 dxdt < / / Vip ■ (X7r/^)dxdt + / r/iptipdx 
o Jnno' Jo Jn Jn 



Therefore using Leibniz's rule and Cauchy-Schwarz's inequality we obtain that 

[ T [ \*p t \ 2 dxdt < d f T [ (|V^| 2 + \^\ 2 )dxdt + C 2 E(0), 
Jo Jnno 1 Jo Jnno 

for some positive constants Ci,C 2 independent of T . Inserting this estimate in (|5.7p we 
get 



2T£(0) <{C + C 2 )E(0) + C 4 [ [ {\V4>\ 2 + \^\ 2 )dxdt, 

Jo Jnno 



(5.8) 



where C4 = max{l,Ci}. To eliminate the last term of this right-hand side we use a 
compacteness/uniqueness argument. Namely assume that (|5.3p does not hold. Then 
there exists a sequence of (ipi)gpf$ solution of (|5.5p with initial data ip® and tpj such 
that 

||V^°|| + ||^|| = 1, (5.9) 



and 



T 



1 



a(x) \Vip t {x,t)\ A dxdt = - (5.10) 



'0 Jn 

Since the system (|5.5p is conservative, the sequence (tpg) is bounded in H 1 (Qt), where 
Qt = 0, x (0, T)). Hence it converges strongly in L 2 (Qt) to ip a weak solution (in 
H 1 (Qt) of (|5.5p with initial data ip°,ip l . But thanks to (|5.8p applied to ^ — tyf, we 
deduce that the sequence (V^gJ^ i s a Cauchy sequence in ^(Qt)- Hence ipi converges 
in H 1 (Qt) to ip, weak solution of (|5.5p and satisfying 



and 

Accordingly 
or quivalently 



IIV/H + 11^11=1, (5.11) 



T r 

a(x) \Vip(x,t)\ 2 dxdt = 0. (5.12) 

Jn 



Vip = in u x (0,T), 
^(z,i) = ip(t) inwx (0,T). 
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But again due to the first identity in (|5,5D . 

^tt(t) = in u x (0,T), 
hence there exists two complex numbers a, 6 such that 

ip(x,t) = at + b in oj X (0,T). 

Now consider the difference 

tjj = if; — at — b. 
Then we see that it a weak solution of (15.51) such that 



i> = inw x (0,T). 
By using Theorem 9.1 of [16] . we deduce that 

= on Qt, 

or equivalently 

^(cc, i) = at + b on Qt- 

But in our situation ip° belongs to and therefore t^(-,t) belongs to H^i^l), for 

all t > 0. As a consequence we deduce that 

= / ip(x,t)dx = (at + b)\n\,\ft > 0. 
Jn 

This imples that a = b = and contradicts (|5.1ip . 
Thus we have shown the following: 

Corollary 5.2. Suppose that 

inf a(x) > 0. 

T/ien the system 15. 1\) is exponentially stable in the sense specified in Theorem \5.1\ 
6 Polynomial stability 

In this section we show that system (jl.ip is polynomially stable if uj contains a neigh- 
borhood of the whole boundary T as in Corollary 15.21 

This result is based on the following result stated in Theorem 2.4 of [5] (see also 
[6j Q3] for weaker variants) . 

Lemma 6.1. A Cq semigroup e tc of contractions on a Hilbert space such that 

p{C) D{ip | € M} = %SL, (6.1) 

satisfies 

\\e tc U \\ <Ct-T\\U \\ v{c) , VU £V(£), Vt > 1, 
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as well as 

\\e tC U \\ < Ct^WUoW^ci), VU eV(£ l ), Vt > 1, 
for some constant C > and /or some positive integer I if and only if 

lim sup — t \\(il3-CY l \\ < oo. (6.2) 

|/3|->oo P 

Lemma 6.2. Assume that co contains a neighborhood of the whole boundary T. 
Then the resolvent of the operator of Ad satisfies condition h6.S\) with I = 3. 

Proof. We use a contradiction argument, i.e., we suppose that (j6.2|) is false with some 
I > 0. Then there exist a sequence of real numbers /3 n —> +oo and a sequence of vectors 
Z n = (u n ,r n ) T in V(Ad) with 

\\ Z n\\ H = 1, (6-3) 

and 

^||(i/3 n -.4d)Z n || w ->0asn-»-oo. (6.4) 
From the definition of .A^ this last property is equivalent to 

f3 l n {if3 n u n + Vr n + au n ) =: f n -> in L 2 (ft) d , (6.5) 

and 

/3^(i/3 n r n + div<u n ) =: 5n -> in L 2 (0). (6.6) 

As usual, taking the inner product of (|6.4j) with Z n , using (|6.3p and the dissipative- 
ness of ^4^) we find 



a \u n \ 2 dx = &((i/3 n - ^d)Z n , Z n ) = o(/3" / ), (6.7) 



where, here and hereafter, a n = o(b n ) means that 



lim t~" = 0- 

n— >oo 0„ 



We multiply (16.5j) by u n , integrate in Q and use (16.31) to find 

(i/3 n \u n \ 2 + Vr n • u n + a|n n | 2 )dx = o(/3~ l ). 
Using (|6.7p and integrating by parts we obtain 

(if3 n \u n \ 2 - r n divu n ) dx = o((3~ l ). 

i 

Using (|6.6p (and (|6.3|) ) we get 

ipn I (\u n \ 2 ~ \r n ?) dx = o(f3~ l ). 
Jo. 
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This shows that 

Nn|| 2 -||r„|| 2 = (/3- (m) ). (6.8) 



The identity (|6.5p implies that 

Vr n = p- l f n -(i/3 n + a)u n , (6.9) 

and therefore 

||Vr n || = 0((3 n ), (6.10) 

where, here and hereafter, a n = 0(b n ) means that there exists C > independent of n 
such that a n < Cb n for n large enough. 

By (|6.9p . we also have 

a\Vr n \ 2 dx<2 / a\i/3 n + a\ 2 \u n \ 2 dx + 2j3~ 21 / a\f n \ 2 dx. 
Jn Jn 



Hence, for I > 2, by using (|6,7p we find that 

a\Vr n \ 2 dx = o{J3 2 n - 1 ). (6.11) 



Now we notice that (|6.5|) and (|6.7|) imply that 

if3 n u n + Vr n = h n := /3~ z .fn - au n , (6.12) 

with 

IIM =o{fi- l/2 ). (6.13) 
We use the Helmholtz decomposition 

h n = V(p n + Xn, (6.14) 

l|Vy n ||n + llXnlln = \\h n \\ 2 

to deduce that 

ll^n|| 2 ,n + IIXn||n = o(^)- (6-15) 
The identities (j6TT2|) and fljQgD yield 

i/3 n «n + Vr n = V<£„ + Xn, 

or equivalently 

i/3 n « n = - Vs n + Xn, (6.16) 

with 

s n = r n -ip n . (6.17) 

As u n is solenoidal, s n is harmonic, and, furthermore since u n ■ n = x n ■ n = onT, 
we automatically have 

Vs n • n = on T. (6.18) 
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Next, by (|6.3p and (|6.15p . one has 

||a„||=0(l). (6.19) 

while (fBTTUD and (pTL5D lead to 

||Vs„|| = 0{p n ). (6.20) 
Now by (|6.6p . ()6.16|) and (|6. 18[) we see that s n satisfies 

ifinSn ~ T^-divVs n = g n (6.21) 

where 

9n ■= Pn l 9n ~ ifin^n- 

Due to (|6.6p and (|6.15p we see that 

ll5n||=o(/?^ /2 ). (6.22) 

At this stage we use the multiplier method to find some properties on s n . First we 
multiply (|6.2ip by — is n and integrate in ft to obtain due to (|6.19p and (|6.22|) 

/ WnSn + "^-div Vs n )s n dx = o{^ 2 ). 

Jn Pn 

By an integration by part and since the boundary term is zero due to (|6.18p . we obtain 

/ (Pn\Sn\ 2 ~ 7^|V S „| 2 ) dx = O^ 2 ). (6.23) 

Jn Pn 

Secondly we take a (smooth) real- valued multiplier m G C 2 (£l) d fixed later on but 
such that m = on V and multiply (|6.21|) by — im ■ Vs n and integrate in to obtain 
due to dOQD and ([6T22]) 

/ {Pns n + ^-div Vs n )m • Vs n dx = o((3l~ 1 / 2 ). (6.24) 
Jn Pn 

Now, in a standard way by means of Green's formula, the first term 

I\ := f s n m ■ Vs n dx 

is transformed into 



n 



h = ~ d k (s n m k )s n dx + / m ■ n\s n \ da 
Jn Jr 

= — Ji — / div m|s n | 2 dx, 
Jn 

recalling that m = 0on the boundary. Hence we have 

28Wi = - / dwm\s n \ 2 dx. (6.25) 
Jn 
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In the same manner we have 



I2 '■= / div Vs„m • Vs n dx = — dkS n dk{m ■ V s n ) dx (6.26) 
Jn Jn 

d k s n d k mjdjS n dx - I 3 , 



where 



n 



h= / d k s n mjd k djs n dx. 
In 



Again an integration by part leads to 



dj(d k s n mj)d k s n dx 
I3 — div m|Vs n | 2 dx. 



9. 



Jn 

Consequently 

= — / div m|Vs n | 2 dx, 
Jn 

and using this expression in (|6.26p we find 



2m 2 = -2K / dksndkmjdjSn) dx (6.27) 
Jn 



+ / div m|Vs n | 2 dx. 
>n 



Taking the real part of (|6.24p and using (I6.25P and (|6.27p we obtain 



-j3 n I divm|s n | 2 dx + — / divm|Vs„| 2 dx 

Pn 



n 



~Sl [ d k s n d k m j dfs n dx = o{p 2 - 1 / 2 ). (6.28) 
Pn Jn 

From (|6.1ip and (j6. 15[) . we see that 

/ a\Vs n \ 2 dx = o(/3 2 " z ), (6.29) 
Jn 

hence for I > 2, Vs n tends to zero on co. This means that we are mainly interested in 
the behaviour of s n outside oj. Therefore we take for m a function with a support far 
from the boundary, namely we take 

771(2) = r](x)x, 
where 77 is a smooth cut-off function such that 

T] = 1 on uj c and r] = in O', 
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where O' is a neigbourhood of T. Since 

dkirij = 5 jk r] + XjdkTj, 

(|6.28p becomes 

-P n [ (di 1 + x-VT])\s n \ 2 dx + ^- [ ((d-2)rj + x- Vr?)|Vs n | 
Jn Pn Jn 

-—3? / dkSndk-qxjdjSndx = o(f3l~ l/2 ). 
Pn Jn 

Assume that the next estimate 

^fj Sn \ 2 d X y =o(^) 

holds. Then combining this estimate with (|6.29p in (|6.30p we obtain 

-Pn* [ r,\s n \ 2 dx + ^-(d-2) [ V \ Vs n \ 2 dx = o(/3^ /2 ). 
Jn Pn Jn 

Coming back to (|6.23p and writting 1 = rj + (1 — rj), we get 

[Wnri\sn\ 2 -^v\Vs n \ 2 )dx = oiPi- 1 ' 2 ) 
Jn Pn 

- /" (^.(1 - 77)|^| 2 - -i-(l - 77)|V^ 
Jn Pn 

Hence (f!T3Tj) and (ET29]) lead to 

/ (M\Sn\ 2 ~ ^~V\ VS„| 2 ) dx = 0{P 1 - 1 ' 2 ) + o($t l ). 

Jn Pn 
This estimate in (I6.32P yields 

r ] \s n \ 2 dx = o((3 1 -^ 2 ) + o((3 2 - 1 ). 

By choosing I > 2 and using again (|6.3ip . we arrive at 

Sn \ 2 dx = o{p l ~ 1 ' 2 ). 

n 

Taking into account (|6.23[) we get 

\Vs n \ 2 dx = p 2 n [ \s n \ 2 dx + o{p 2 n - l l 2 ), 



and (I6.34p finally leads to 

' \Vs n \ 2 dx = o(f3t l/2 )- 
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Coming back to our original variables, for r n we have by ()6.34j) and (|6,15p that 

IMI 2 = o{&- l i% 

while for u n using (|6. 16|) . (|6.35p and (|6.15p 

\\u n f = o(^). 

In conclusion for I > 3, ||r n || and \\u n \\ tend to zero which contradicts (|6.3|) . 

It remains to prove (|6.3ip . For that purpose, we first show that the mean of s n on 
T tends to zero. Fix a function h £ (C 2 (Sl)) d such that supp h <Z uj and 

ft ■ n = 1 on T. 

Such a function exists by Lemma 1.3.1 of [13] since we assume that Q, has a C 3 -boundary. 
We then multiply (|6.5p by h and integrate in f2 to get 



((i/3 n + a)-u n + Vr n )ft dx = o(/3 n l ). 

i 

Hence integrating by parts, we obtain 

((if3 n + a)u n h - r n divh) dx + / r n da = o((3~ l ). (6.36) 



We still need to transform the term 



r n div hdx. 

Jn 

By using (I6.6p we may write 

/ r n divhdx = — i \ (/3~^ l+1 ^ g n — div u n )div hdx. 
Jn Jn 

Again by Green's formula and the fact that u n ■ n = on T we obtain 

/ r n div hdx = —i / (/3~^ +1 ^ n div h + u n - Vdiv h) dx. 
Jn Jn 

Using this identity in (|6.36p we obtain 

r n da = — \ {{i@ n + a)u n h — if3~^ +l ^ g n d\v h — iu n ■ Vdiv/i) dx + o(f3~ l ). 



Using (|6.7p we get 



j^r n do = o{p l n - l l 2 ). (6.37) 
But a standard Poincare type inequality implies that 

(J |r n | 2 ^ 2 <C{(J \Vr n \ 2 d x y + f^r n d<j , 
for some C > depending only on us. Hence using (|6.37p and (|6.1ip yields 



\r n \ 2 dxj =o(f3 1 n - 1 ' 2 ). (6.38) 
This estimate and (|6.15p leads to (|6.3ip . □ 
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This Lemma and Lemmas 13.31 and 13.41 show that the restriction of the operator Ad 
to Hq satisfies the hypotheses of Lemma 16.11 with / = 3. Therefore we have obtained 
the next main result. 

Theorem 6.3. Assume that oo contains a neighborhood of the whole boundary V then 
the system (jl.ip with an initial datum (u°,r°) in Hq n D(Ad), Hq = Ker[„4,d]~S is 
polynomially stable, namely there exists C > such that 

E(t)<ct- 2 / 3 \\(v°y)\\l iAd) yt>o. 
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